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IMPROVEMENTS RELATING TO VELOCITY EXTRACTION 

This invention concerns improvements relating to velocity extraction in 
the field of radiation pulse echo detection. More particularly but not exclusively, 
this invention concerns extraction of the velocity of a target from the returns. 

Radiation pulse echo detection systems, such as radar, transmit a set of 
one or more coherent strings of pulses (coherent bursts) which are reflected by 
objects. The echoes of pulses are used to detect and locate distant objects. 

Conventionally, a surveillance radar would estimate the target radial 
velocity using multiple estimates of the target range taken at different times. It 
is now common practice to filter the returns to remove clutter returns (those 
returns from items not of interest), leaving any returns from moving items of 
interest relatively unaffected. This is known as Moving Target Detection (MTD) 
or Moving Target Indication (MTI). Such schemes can be modified using 
multiple filters to obtain some measure of target radial velocity but suffer from 
problems due to large clutter returns 'spilling 1 into adjacent filters, thereby 
producing erroneous velocity measurements. 

One of the main tasks of modern radars and sonars, is to identify and 
track moving targets. The accuracy of tracking is greatly enhanced if the range 
ambiguity and radial velocity of the target input plots are known. Only plots with 
matching range ambiguity and velocity will then be associated with those from 
previous measurements thereby significantly reducing the probability of mis- 
association, track seduction and false track rate. 

Whilst traditional filtering methods remove the clutter from the in-phase 
(I) and quadrature (Q) components (the first being in phase with the transmitted 
signal and the second in quadrature with the transmitted signal) and typically 
return the target amplitude, they do not easily output the target radial velocity 
nor the range ambiguity of the target. 

The present invention aims to overcome or at least substantially reduce 
some of the above mentioned problems. 

According to a first aspect of the invention, there is provided a method of 
extracting the target radial velocity from one or more coherent radiation pulse 
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bursts comprising the steps of: (a) receiving radiation echo returns from a 
scene; (b) processing the returns into in-phase (I) and quadrature (Q) 
components; (c) measuring returns at intervals to provide sampled data; (d) 
applying a function to the l-Q returns; (e) modifying the function to match the 
sampled data as a function of velocity in a predetermined fashion; and (f) 
determining the target radial velocity in dependence upon said modification step 
of the function. This method advantageously produces accurate velocity 
measurements and is a real time process, fast enough to work at, for example, 
radar type timescales with data rates of a few MHz. The method is applicable 
to various radiation pulse echo detection systems that use coherent bursts for 
clutter rejection purposes. 

Preferably, step (d) comprises fitting a curve to the l-Q returns (using 
known system parameters) and step (e) comprises optimising the fit to the 
sampled data as a function of velocity in a least squares fashion to give the 
target radial velocity. 

Other target parameters, such as target amplitude, range ambiguity and 
azimuth, may then be derived from the fit to the sampled data. 

Conveniently, the model of clutter return, used in the fitting of a curve to 
the l-Q returns, is a low order polynomial in I and Q. This is a simple model that 
provides a good fit to clutter return. Alternatively, if desired, the model of clutter 
return could be a different analytic function. 

The returns may be measured at non-equi-spaced intervals. The use of 
non-equi-spaced intervals provides unambiguous velocity measurement 
whereas sampling at a regular rate can lead to multiple solutions. 

The transmission frequency may be altered between bursts. This also 
gives unambiguous velocity measurement. 

Each burst preferably consists of multiple pulses transmitted at irregular 
intervals (non-constant pulse repetition interval bursts). Digital Fourier 
transform techniques previously used for the same purpose of target parameter 
extraction cannot easily be used with such in-burst stagger. 
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The bursts are optionally internally coherent but mutually incoherent. 
This gives an improved measure of target velocity, amplitude, range ambiguity 
and azimuth. 

The fit residues from adjacent range cells in which a target was detected 
5 may be summed to obtain the target velocity and other parameters. This 
increases the probability of detecting targets split between range cells. 

Conventional MTI/MTD filtering may be carried out before applying a 
function to the l-Q returns in which a target was detected. Post-detection 
processing of the data requires a lower processing capacity than pre-detection 
10 processing. 

Optionally, post-detection processing is used in combination with pre- 
detection processing to limit the range cells processed to increase efficiency. 

It is to be appreciated that the present invention may be embodied in 
software. Accordingly, the present invention extends to a computer program 
15 element comprising program code for configuring a programmable device 
apparatus or system to implement the above described method. Suitably, the 
computer program is stored on a carrier medium. 

Further, the present invention extends to a data processing system or 
apparatus adapted and arranged to implement the above described method. 

20 Preferably, there is provided a data processing system comprising a 

transmitter; an antenna; a receiver; signal processing means; an I and Q 
component splitter; an analog-to-digital converter, and processing means to fit a 
predetermined function to the I and Q components. 

The invention will now be described by way of example and with 
25 reference to the accompanying drawings, in which: 

Figure 1 is a flowchart illustrating the steps involved in a method 
embodying the invention. 

Figure 2 shows a typical return in l-Q space. 

Figure 3 is a flowchart of a variant of the method shown in Figure 1 . 



i 
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Figure 4 shows an alternative embodiment of the invention. 

Figure 5 is a flowchart illustrating a variant of the method shown in Figure 

4. 

Figure 6 is a flowchart of an alternative embodiment of the invention. 

Figure 7 shows a variant of the method shown in Figure 6. 

Figure 8 shows the optional steps for extraction of range ambiguity. 

Figure 9 is a schematic view of a data processing system in accordance 
with one aspect of the invention. 

Figure 1 shows the steps of a method embodying the invention. As 
shown, a coherent burst is emitted (step 102) from a radar antenna and 
reflected by objects in its path. The returns from a coherent burst are split (step 
104) into two mutually orthogonal components, the in-phase (I) and quadrature 
(Q) components. Measurements are made (step 106) of the return signal 
(returns) at predetermined time intervals. The returns from the scene constitute 
noise (clutter returns) and target returns, if a target is present. 

The returns from a target will typically constitute a helix whose long axis 
is centred along the time axis. The clutter returns cause the centre of the helix 
to lie along the clutter trajectory. For stationary clutter this is a straight line 
parallel to the time axis; for moving clutter the line is curved, but substantially 
parallel to the time axis. 

It is therefore assumed that in l-Q space the return signal consists of a 
noise-contaminated helix. Figure 2 shows an exemplary return signal 24 in l-Q 
space. The clutter trajectory 26 is substantially parallel to the time axis, 
although curved due to moving clutter. A target helix 28 is imposed onto the 
clutter trajectory. The motion of the target imposes a changing phase shift on 
the signal which is used to estimate the radial velocity from the rotation rate of 
the sample vectors, i.e. the helix pitch is proportional to the target radial 
velocity. The radius of the helix is proportional to the target cross-section at the 
transmitted frequency, and the helix axis lies on a clutter return curve that can 
be represented by a low order polynomial in time. 
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Returning to Figure 1, a curve constituting a simple low order polynomial 
in I and Q model of the clutter return, and a helical model of the target return is 
therefore fitted (step 108) to the sampled data to describe the returns assuming 
a target with a guess velocity is present together with clutter. 

The fit to the sampled data is optimised (steps 108a - e) in a least 
squares fashion to minimise the error value, or residue, e, given by 
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= A b is the best estimate of the mean amplitude of the target at time 



/, seen in burst b ; 

Isample b mi Qsample b m are the in-phase and quadrature components of the 
measured signal for pulse m in burst b ; 

h t n ■ Qb, n are the best estimate of the components of the measured signals due 
to the clutter seen in burst b , and 
0 6 is the phase at the start of burst b . 

The best fit target radial velocity is obtained (step 110) and other target 
parameters are derived (step 112), for example the target amplitude is extracted 
from the radius of the helix. The return strength is compared (step 114) with a 
threshold target return strength which is chosen to produce a desired probability 
of false detection. 

Integrating several bursts together improves the detection probability for 
a target. The number of bursts over which summing occurs is referred to as the 
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detection or integration window. For rotating radars where target returns are 
smeared by the rotation of the antenna, the integration window is usually set to 
the number of bursts transmitted in the time the radar sweeps a beamwidth. 

Since targets may exceed one range cell in size or bridge range cell 
boundaries the accuracy of the velocity extraction process can be increased if 
any adjacent cells that also contain returns from the target are included along 
with those cells in the normal integration window. Similarly, since the beam 
shape of the radar is smeared in azimuth, the same target will be detected at 
the same range by several bursts. These detections in adjacent bursts at a 
common range are also grouped together as belonging to the same plot. This 
clustering of adjacent detections is known as 'plot extraction'. The various 
methods to do this clustering, whether in range then azimuth or in azimuth then 
range, and the methods for grouping detections at common ranges and 
azimuths but differing, though adjacent, elevations for three-dimensional radars, 
are known as 'plot collapsing' (step 116). These methods enable a single plot 
to be formed, corresponding to the best estimate of the target position, rather 
than a cluster of detections (partial plots). The most common methods group in 
range, providing a plot start and stop range in each burst, then group any 
azimuthally adjacent overlapping detection runs into the cluster. The plot 
position is then set to the centre of the cluster. The centre position can be 
found in a number of different ways, for example, the 'centre of mass 1 of the 
cluster or from halving the position of the extreme position of the partial plots 
forming the plot in range and azimuth. Alternatively, the partial plots are passed 
to a track extractor for outputting parameters. Track extractors make use of 
information from previous scans, knowledge of target range or velocity 
ambiguities to improve the accuracy of the information extraction and to resolve 
any plot track association conflicts (when two or more targets are too close 
together to be resolved by conventional plot collapsing methods but where the 
knowledge that two or more targets are expected close together with, for 
example, different velocities enables the returns to be separated into two more 
more clusters separated in velocity space). 



XA1595WOspecdoc 



-7- 

For each partial plot (i.e. detection) in the cluster, the ratio of the mean 
residue to the minimum residue is the mean return power-to-noise ratio for the 
target in that cell. By weighting the radial velocity from each cell in the cluster 
by the mean return power-to-noise in that cell, summing all the weighted 
velocities and dividing by the sum of the weights, the weighted mean radial 
velocity is obtained (step 118). This is a better estimate of the radial velocity 
than a simple mean of the velocities as it gives greater weight to those 
velocities that have better signal-to-noise ratios and thus accuracies. 

The burst parameters are averaged within the clusters to give the cluster 
parameters (step 120) and the parameters are outputted (step 122). 

Optionally, the residues from adjacent range cells with bursts that cross 
the detection threshold are summed (see Figure 3; where the final two digits of 
the numbered steps are the same as those in numbered steps beginning with a 
T in Figure 1 the steps are equivalent, but not necessarily identical). As each 
burst is processed, detections are either declared or not for each range cell and 
the residues stored for those bursts with detections (step 330). The residues for 
each velocity step are summed (step 332) and the minimum residue found (step 
334). 

The target radial velocity and other parameters are obtained (steps 336, 
338) from the new summed residue, i.e. summation takes place across the 
partial plots, not across the integration window. The original (partial plot) 
detection threshold for target strength at step 314 can then be quite low with a 
high plot detection threshold applied (step 340) to the summed data to provide 
the required probability of false plot detection. The lower original threshold at 
step 314 enables smaller targets to be detected, although with a higher 
probability of false alarm. The higher post interaction threshold at step 340 re- 
establishes the lower probability of false detection. As more bursts are 
integrated together, true targets that are split between range cells, for instance, 
integrate up and have a higher probability of detection. 
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ln Figures 3, 4, 5, 6 and 7, where the final two digits of the numbered 
steps are the same as those in numbered steps beginning with a '1' in Figure 1 
the steps are equivalent, but not necessarily identical. 

In a further embodiment of the invention, the integration process is 
carried out post-detection. While having a lower sensitivity, post-detection 
processing of the data requires a lower processing capacity than pre-detection 
processing. In the method shown in Figure 4, conventional MTI/MTD filtering is 
carried out (step 444) on the I and Q components of the echo returns. In both 
MTI and MTD filtering, the filters are equivalent to multiplying the returns within 
a burst by a set of weights and summing the results, the weights being chosen 
to reject the expected clutter, whilst allowing the returns from moving clutter to 
sum constructively. The I and Q components for the bursts in each range cell 
where detections occur are stored (step 446) before conventional plot 
collapsing methods are used (step 416) to detect clusters. I and Q dependent 
terms of residue are calculated (step 408a) for each burst in a group of adjacent 
partial plots (in each 'cluster 1 of detections) as a function of velocity. For each 
burst the residue is calculated (step 408b) and stored (step 408c) for each 
velocity step. The residues are then summed (step 408d) across all bursts 
included in the plot. The minimum residue is found (step 408e) and the target 
radial velocity and other parameters determined (steps 410, 412) before being 
outputted (step 422). 

Alternatively, the plot collapsing (step 516) may be carried out after the 
velocity extraction processing (see Figure 5). In this case, each partial plot 
contributes a velocity measure, the weighted mean radial velocity is then 
calculated (step 518) and the parameters outputted (step 522). This process is 
more computationally efficient (the velocity being extracted from only those 
range cells where detections occur) but less accurate than the process 
illustrated in Figure 4 where the velocity is extracted from data for all range cells 
in the cluster that forms the plot. 

Clearer identification of minima is obtained by enhancing the difference 
between minima and peak values. This can be done by summing together all 
curves but this requires a great amount of computing power. The use of post- 
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detection processing reduces the computing power required. The post- 
detection processing method involving plot collapsing before the calculation of 
residues (Figure 4) is preferred to the post-detection processing method 
involving the calculation of the weighted mean radial velocity (Figure 5) as more 
bursts are included in the calculation so the likelihood of target detection is 
increased. 

It is more efficient to discard those cells which clearly only contain noise, 
before carrying out the velocity extraction processing. This extra level of 
processing can reduce the processing load further downstream by significantly 
reducing the number of cells which receive the full pre-detection processing. 

By applying a simple MTI filter with a low threshold (as compared to the 
threshold used in the processes illustrated in Figures 4 and 5) to the summed 
returns in all cells, the large number of cells where no target is illuminated and 
that clearly only contain clutter or noise and that have a very low probability of 
containing a target as seen from their spectral content are eliminated. In this 
manner, more real targets are detected and less false targets are found, 
compared to the processing methods shown in and described with reference to 
Figures 1,3,4 and 5. Examples of such methods are illustrated in Figures 6 
and 7. 

In the method shown in Figure 6, conventional MTI/MTD filtering is 
carried out (step 650) using a reduced threshold, on the I and Q components of 
the echo returns. The use of a reduced threshold increases the likelihood of 
target detection. The I and Q components for the range cells with bursts where 
detections occur are stored (step 646) before conventional plot collapsing 
methods are used (step 616) to detect clusters. I and Q dependent terms of 
residue are calculated (step 608a) for each burst in a group of adjacent partial 
plots (in each 'cluster 1 of detections) as a function of velocity. For each burst 
the residue is calculated (step 608b) and stored (step 608c) for each velocity 
step. The residues are then summed (step 608d) across all bursts included in 
the plot. The minimum residue is found (step 608e) and the target radial 
velocity and other parameters determined (steps 610, 612). The target return 
strength is compared to a threshold value (step 614) and the parameters are 
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outputted if the target return strength is greater than the detection threshold 
(step 652). The threshold value is chosen to maintain the false detection rate at 
the required level, previously increased by the use of a reduced threshold 
during filtering (step 650). 

5 In an alternative method, the plot collapsing (step 716) may be carried 

out after the velocity extraction processing (see Figure 7) and the burst 
parameters are averaged within clusters to give the cluster parameters (step 
754). The target return strength is compared to a threshold value (step 714), 
the mean weighted radial velocity calculated (step 718) and the parameters 
10 outputted if the target return strength is greater than the detection threshold 
(step 752). The process shown in Figure 7 is more computationally efficient but 
less accurate than the Figure 6 process. 

When a coherent burst is transmitted objects will be illuminated at ranges 
that are sufficiently far away that the echoes from any particular pulse are not 
15 received back at the radar until after one or more further pulses have been 
transmitted. These longer-range returns are ambiguous in range; the radar will 
normally interpret them as coming from a close range as calculated from the 
time interval from the time of the reception and the transmission of the last 
pulse prior to reception. The problem can often be avoided by using in-burst 

20 stagger (i.e. the pulse repetition interval (PRI), the time between transmitted 
pulses, is not constant); the echoes from ambiguous range targets are spread 
over a number of different ranges and often do not integrate up to a sufficient 
level to give a high probability of detection. For constant PRI bursts the 
ambiguous returns from any one target will, within that burst, all lie at the same 

25 erroneous range and will coherently sum to give a larger return, which may then 
be detected. For multiple bursts at different frequencies but with common PRI 
these ambiguous range targets then integrate across the bursts to give false 
detections. These false detections, i.e. the reporting of a real target at an 
incorrect location, are problematic for many MTI and MTD radars. The problem 

30 is exacerbated when multiple bursts with different PRIs are integrated together. 
In this case, whilst small ambiguous range targets do not integrate across the 
bursts as they appear at different locations in each burst, larger targets with 
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sufficient return strength for detection to occur on the basis of a single burst 
then appear as detections in up to as many different locations as there are 
bursts integrated. Current methods of dealing with this problem, commonly 
known as ambiguous range trap (ART) methods suffer from a number of 
problems, most notably the blanking of range cells where ambiguous targets 
appear to be and therefore the erroneous deletion of cells with real targets 
present. 

The method in accordance with this invention provides a means, not only 
of reliably identifying when the returns in the burst come from an ambiguous 
range target, but also of obtaining the order of ambiguity of the target. 

The target range ambiguity can be estimated by considering whether a 
better fit to the data would be obtained by assuming that one or more of the 
initial returns do not lie on the helix. A target is a zero order ambiguous target 
(i.e. an unambiguous range target) if the first transmitted pulse is received 
before the second pulse is transmitted. Thus, for a zero order ambiguous 
target, the return from the first pulse would contain a return from the target and 
lie on the helix. An ambiguous range target (one with an ambiguity order higher 
than zero) would not have any signal from the target contained in the first 
received pulse and thus the received signal would lie on the axis of the helix. 
The ambiguity order n is given by the number of first received pulse signals 
lying on the axis of the helix, only the subsequently received signals lying on the 
helix itself. 

When it is known that only a single or at least a very sparse set of targets 
is present, or where the bursts being integrated together all have the same PRI 
(although they may have different transmission frequencies) it is possible to 
carry out a simple minimising of the residue with respect to ambiguity and 
velocity. 

In the latter case, when all the bursts use a common PRI, the previously 
described processing in any of the variants may be used with the addition that 
the residues are calculated for each of the possible range ambiguities of the 
target. The target velocity and range ambiguity is found simply by determining 
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the velocity and range ambiguity for which the residue is a minimum. In the 
former case, where the PRIs differ from burst to burst but it is known that only 
one, or at least very few, targets are present, the process is slightly more 
complicated in that the returns from each burst are first replicated to all possible 
range ambiguities. The unfolded returns from different bursts are then 
realigned over the possible ambiguous ranges so the residues from returns at 
the ambiguous ranges in all bursts can be calculated and summed together for 
each possible ambiguity. At the correct range the returns from all the bursts 
align, and therefore the minimum will be deepest. Summation of residues at all 
possible ambiguous ranges will give the deepest minimum thereby solving the 
range ambiguity and target velocity. Having determined velocity and range 
ambiguity the processing then follows that of the previously described methods. 

For the more normal case where the PRIs differ from burst to burst and 
the potential number of targets is not very small, a more convoluted process is 
required to extract the range ambiguity. The simplest approach to extracting 
the range ambiguity is to extract it on a burst-by-burst basis. The minimum 
residue for each range cell, or at least those where there is a possibility that a 
target is present, with respect to both velocity and range ambiguity is 
determined. The range ambiguity at which this minimum occurs is taken to be 
the range ambiguity order of the return. From the apparent range and the order 
of the ambiguity the actual range of the return can be determined. 

With ambiguous range returns, a choice of methods equivalent to the 
various methods of extracting the radial velocity from an unambiguous return is 
available. The returns from the bursts within an integration window at a 
common (ambiguous) range can be processed to give a maximum likelihood 
velocity and target size, or any other of the available parameters extracted 
either by a simple weighted mean approach or by re-computing the residues for 
the whole set of returns. For long coherent bursts this approach suffers from 
the difficulty that the difference in depth of the residue between adjacent orders 
of ambiguity is small and therefore it is possible that for relatively weak targets 
the wrong order of ambiguity may initially be chosen. Errors of this sort are 
typically only of one or two ambiguities; by unfolding the target to just those 
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ambiguities around the most likely ambiguity order and summing the residues 
across the bursts in the integration or plot window (depending on the method of 
combining the bursts used), then seeking the deepest minimum from this limited 
set of possible target locations, the most likely velocity, position and other target 
5 parameters can be extracted from the data. The most likely azimuth may 
possibly also be extracted. 

For those processes that involve pre-detection processing, such as 
illustrated in Figures 1 and 3, the range ambiguity of each burst is determined 
(step 112 and 312 respectively) and bursts with returns showing targets that 

10 form clusters in range and/or range rate space, as determined from the 
extracted range ambiguity, are then combined. In the former case, where 
clusters are formed in range, the velocity is determined, on the assumption that 
only one target is present, by processing all bursts with returns above the 
detection threshold at that range to extract the best fit velocity and mean return 

15 strength. In the latter case, where clusters are formed in range rate space, only 
bursts with velocities that unfold to a common velocity are included in the 
processing. This works best when there is a high probability of two targets with 
differing velocities, spatially unresolved, at the same range. 

For processes that involve post-detection processing, the bursts included 
20 are those involved in the detection, irrespective of whether their range 
ambiguities match, i.e. the range ambiguity is solved at the same time as the 
velocity is extracted from all contributing bursts. In the processes shown in 
Figures 4 to 7, steps 408b to e, 508b to e, 608b to e and 708b to e are replaced 
with steps 8b 1 to 8e\ shown in Figure 8. For each burst and possible range 
25 ambiguity, the residue is calculated (step 8b') for each velocity step. The 
residues are stored (step 8c') and summed (step 8d') (across bursts included in 
the plot) for the possible ambiguities for each detection. The minimum residue 
is then found (step 8e') across both range and velocity (at each range ambiguity 
the velocity corresponding to the minimum residue will be different). 

30 Alternatively, in Figures 4 to 7, residues can be determined in all 

range/azimuth cells out to the maximum possible ambiguity for all range cells 
with detections for all bursts over some sector. The use of a fixed sector width 
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overcomes the problem that a cluster at one ambiguity may partially overlap 
with a cluster at another, preventing automatic decisions on processing ranges. 
A model is then set up with all possible target position/velocity combinations 
over the sector space and the residues corresponding to those combinations 
5 found. The deepest minimum over the multidimensional (2n) space (positions 
and velocities for each of the n possible targets) is selected. The velocity and 
range corresponding to the minimum residue is then determined. 

When extracting range ambiguity, two types of false result can occur. A 
false velocity may be outputted and false targets may be 'detected' due to 
10 ghosts (where returns from two or more targets at different ranges seen in 
bursts with different pulse repetition intervals unfold onto the same range with 
differing ambiguities resulting in an erroneous detection). Both false outputs 
can be reduced by including as many bursts where the target is present as 
possible since true targets integrate up while false targets do not. 

When range ambiguity is extracted using pre-detection processing, the 
number of cells in which a false alarm can occur is increased, due to the 
unfolding of ranges to the full ambiguous range limit. Therefore a higher 
threshold is needed to maintain the required false detection rate and reduce the 
possibility of the residue minima being due to ghosts. This higher threshold is 
applied after plot collapsing to maintain the maximum sensitivity. The minima in 
the set of cells corresponding to potential ghosts are compared, the deepest 
minimum being the most likely solution to the target's radial velocity and range 
ambiguity. 

If a target is range ambiguous and the range ambiguity is not extracted, 
25 the fit to the data is not as accurate as possible, giving rise to an inaccurate 
velocity extraction. 

Tabulation of the error values with respect to velocity and/or range 
ambiguity gives a measure of confidence in the relevant measurement. The 
sharper the minimum, the narrower the uncertainty in the measurement. The 
30 deeper the minimum with respect to the next deepest, the more confidence 
there is in the result. 



20 
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The azimuth of the target (the angular direction of the target from a 
defined direction, e.g. north, in a horizontal plane) is then calculated (for rotating 
radars only) by taking the ratio of the rate of change of target returns strength to 
returns strength for each coherent burst included in the process and fitting the 
5 resultant normalised gradients to the beam shape. The resultant ratios are 
used to derive the polynomial 



the solution for which gives the time at which the antenna pointed at the target. 
From that time, the azimuth of the target is determined as the azimuthal 
10 direction of the antenna at time / . 

The crossing rate of the target, i.e. the rate at which the target crosses 
the field of view of the antenna, can also be determined from the beam shape. 
This can be combined with the radial velocity to give the target velocity. 

An embodiment in accordance with another aspect of the invention is 
15 shown in Figure 9. A data processing system 60 comprises a transmitter 62; an 
antenna 64; a receiver 66; signal processing means 68; an I and Q component 
splitter 70; an analog-to-digital converter 72; a digital filter 74 and processing 
means 76 to fit a predetermined function to the I and Q components. The 
transmitter 62, connected to processing means 76 (connection not shown) or 

20 other processing means (not shown), emits coherent radiation bursts from the 
antenna 64. Radiation echo returns of the pulse bursts are received from a 
remote scene by receiver 66. Signal processing means 68 may include, for 
example, RF filtering to remove out of band signals, automatic gain control 
(AGC) to reduce the gain when strong signals are received to prevent overload, 

25 beam forming to focus the radar to look in one specific direction, jamming null 
steering to reduce the antenna gain in the direction of jammers or pulse 
compression to allow the transmission of long low mean power signals that can 
be compressed on reception into a short high power return. The echo returns 
are processed into I and Q components by the splitter 70 either before or after 

30 processing by the analog-to-digital converter 72. The returns are measured at 
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intervals to provide sampled data and conventional MTI/MTD filtering is applied 
by the digital filter 74. Preferably, the digital filter 74 comprises a digital signal 
processor (DSP) or conventional high speed microprocessor. The filter applied 
can be, for example, a Fast Fourier Transform (FFT) or weighted sum. 
Processing means 76 then applies a predetermined function to the l-Q returns 
and modifies the function to match the sampled data as a function of velocity in 
a predetermined fashion. The target radial velocity is determined in 
dependence upon the modification step of the function. 

Having now described various embodiments of the invention, numerous 
modifications will become apparent to the skilled person. For instance, 
sampling returns at a regular rate can lead to multiple velocity solutions. It is 
therefore advantageous to measure the returns at non-equi-spaced intervals to 
provide unambiguous velocity measurement. This is further aided by the 
alteration of transmission frequency between bursts. It is preferable for the or 
each burst to consist of multiple pulses transmitted at irregular intervals (non- 
constant pulse repetition interval bursts). Digital Fourier transform techniques 
previously used for the same purpose of target parameter extraction cannot be 
used with such in-burst stagger; if filter techniques are required, to reduce 
processing load for example, then MTI processing is required. The use of 
internally coherent but mutually incoherent bursts gives an improved measure 
of target velocity, amplitude, range ambiguity and azimuth. 

To improve the processing time, the velocity independent components of 
the residue may be calculated once, prior to the computation of the residue, 
rather than each time the residue is computed for the different measurements. 
Similarly, the data independent components of the residue can be pre- 
calculated to significantly reduce the processing required to obtain the residue. 

While a polynomial in time / is preferred for describing the returns, other 
analytic forms in orthogonal function space may be used. 

There are various methods of finding the minimum residue. These 
include the computationally intensive method of 'over sampling' in velocity 
space and choosing the velocity corresponding to the smallest value (as 
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described above). A more computationally efficient method calculates the 
residues at the 'Nyquist' frequency (the cut-off frequency above which a signal 
must be sampled in order to be able to reconstruct it) and uses an interpolation 
technique about the minimum calculated point to determine the actual minimum 
and corresponding velocity. 

Where two or more minimum residues are found to be of similar value 
additional processes can be used. This processing may either flag that there is 
a potential problem with the extracted parameters (including velocity) or pass 
the parameters and velocities for all the relevant minima to subsequent radar 
processing (track extraction). This approach minimises the errors due to 
velocity ambiguities and multiple targets in the same range cell on the extracted 
velocity. 



